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Analyses of laser ranges to the Moon provide increasingly stringent limits on any violation of 
the equivalence principle (EP); they also enable several very accurate tests of relativistic gravity. 
These analyses give an EP test of A(Mg/Mi)ep = (-1.0 ±1.4) x 10" 13 . This result yields a strong 
equivalence principle (SEP) test of A(Mg/Mi)sep = ( — 2.0±2.0) x 10~ 13 . Also, the corresponding 
SEP violation parameter n is (4.4 ± 4.5) x 10~ 4 , where r\ = 4/3 — 7 — 3 and both j3 and 7 are 
post-Newtonian parameters. Using the Cassini 7, the rj result yields /3 — 1 = (1.2 ± 1.1) x 10~ 4 . The 
geodetic precession test, expressed as a relative deviation from general relativity, is K gv — — 0.0019± 
0.0064. The search for a time variation in the gravitational constant results in G/G = (4±9) x 10~ 13 
yr" 1 ; consequently there is no evidence for local (~1 AU) scale expansion of the solar system. 

PACS numbers: 04.80.-y; 04.80.Cc; 95.10.Eg; 95.55.Pe; 96.20.Jz 



I. INTRODUCTION 

Einstein's general theory of relativity (GR) began its 
empirical success in 1915 by explaining the anomalous 
perihelion precession of Mercury's orbit, using no ad- 
justable theoretical parameters. Shortly thereafter, Ed- 
dington's 1919 observations of star lines-of-sight during 
a solar eclipse confirmed the doubling of the deflection 
angles predicted by GR as compared to Newtonian and 
Equivalence Principle (EP) arguments. Following these 
beginnings, the general theory of relativity has been ver- 
ified at ever-higher accuracy. Thus, microwave ranging 
to the Viking landers on Mars yielded a ~0.2% accuracy 
via the Shapiro time-delay |l( . Spacecraft and planetary 
radar observations reached an accuracy of ~0.15% [J. 
Lunar Laser Ranging (LLR) has provided verification of 
GR improving the accuracy to ~0.05% via precision mea- 
surements of the lunar orbit 0, 0, 0, El • The astromet- 
ric observations of the deflection of quasar positions with 
respect to the Sun performed with Very-Long Baseline 
fnterferometry (VLB1) improved the accuracy of the GR 
tests to ^0.045% [1,13 • The time-delay experiment with 
the Cassini spacecraft at a solar conjunction has tested 
gravity with a remarkable accuracy of 0.0023% [lfjj in 
measuring the deviation of the parametrized post-Newto- 
nian (PPN) parameter 7 from its GR value of unity. 

The accuracy of the Cassini result opens a new realm 
for tests of gravity in the solar system, especially those 
motivated by the progress in scalar-tensor theories of 
gravity. In particular, the recent work in scalar-tensor 
extensions of gravity that are consistent with present cos- 
mological models [3 13 EI E3 predicts deviations of 
the parameter 7 from the general relativistic value at 
levels of 10~ 5 to 10~ 7 . This prediction motivates new 
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searches for very small deviations of relativistic gravity 
from GR in the solar system and provides a robust the- 
oretical paradigm and constructive guidance for experi- 
ments that would decrease the uncertainty of PPN pa- 
rameters. Thus, in addition to experiments which probe 
the parameter 7, any experiment extending the accuracy 
in measuring parameter [3 is also of great interest. Today 
LLR, the continuing legacy of the Apollo program, is well 
positioned to address this challenge. 

Motivated by the remarkable accuracy of the Cassini 
test for the PPN parameter 7 this Letter reports 
the results of recent LLR analyses using data to April 
2004. The focus is on the improvement of accurate LLR 
gravity experiments, especially the tests of EP, Strong 
Equivalence Principle (SEP), PPN parameter (3, and G. 



II. FUNDAMENTAL PHYSICS WITH LLR 

LLR has a history [4] dating back to the placement of a 
retrorefiector array on the lunar surface by the Apollo 11 
astronauts. Additional reflectors were left by the Apollo 
14 and Apollo 15 astronauts, and two French-built reflec- 
tor arrays were placed on the Moon by the Soviet Luna 
17 and Luna 21 missions. LLR accurately measures the 
time of flight for a laser pulse fired from an observatory 
on the Earth, bounced off of a corner cube retrorefiector 
on the Moon, and returned to the observatory. For a gen- 
eral review of LLR see Dickey et al. . A comprehensive 
paper on LLR tests of gravitational physics is Williams 
et al. Q. A recent test of the EP is in Anderson and 
Williams [|| and other gravitational physics tests are in 
Williams et al. ^5|. An overview of the LLR gravita- 
tional physics tests is given by Nordtvedt ^(|. Reviews 
of various tests of relativity, including the contribution 
by LLR, are given in Will . 

The LLR measurements contribute to a wide range 
of scientific investigations S II 13 EES, and are to- 
day solely responsible for the production of the lunar 
ephemeris. On the fundamental physics front, LLR pro- 
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vides the only current solar system means for testing the 
SEP — the statement that all forms of mass and energy 
contribute equivalent quantities of inertial and gravita- 
tional mass. In addition, LLR is capable of measuring 
the time variation of Newton's gravitational constant, G, 
providing the strongest limit available for the variability 
of this "constant." LLR can also precisely measure the 
de Sitter precession — effectively a spin-orbit coupling af- 
fecting the lunar orbit in the frame co-moving with the 
Earth-Moon system's motion around the Sun. Finally, 
current LLR results are consistent with the GR gravit- 
omagnetic effect on the lunar orbit within 0.1% of the 
predicted level 0,0]. Thus, the lunar orbit is a unique 
laboratory for gravitational physics where each term in 
the relativistic PPN equations of motion is verified to a 
very high accuracy. 



A. Equivalence Principle Tests 

The Equivalence Principle, the exact correspondence of 
gravitational and inertial masses, is a central assumption 
of general relativity and a unique feature of gravitation. 
It is this principle that predicts identical accelerations 
of compositionally different objects in the same gravita- 
tional field, and also allows gravity to be viewed as a 
geometrical property of spacetime-leading to the general 
relativistic interpretation of gravitation. EP tests can 
therefore be viewed in two contexts: tests of the foun- 
dations of the standard model of gravity (i.e. general 
relativity), or as searches for new physics because, as 
emphasized in H 13 , 

almost all extensions to 
the standard model of particle physics generically predict 
new forces that would show up as apparent violations of 
the EP. 

The weak form of the EP (the WEP) states that the 
gravitational properties of strong and electro-weak in- 
teractions obey the EP. In this case the relevant test- 
body differences are their fractional nuclear-binding dif- 
ferences, their neutron-to-proton ratios, their atomic 
charges, etc. GR and other metric theories of gravity 
assume that the WEP is exact. However, extensions 
of the standard model of particle physics that contain 
new macroscopic-range quantum fields predict quantum 
exchange forces that generically violate the WEP be- 
cause they couple to generalized "charges" rather than 
to mass/energy as does gravity 0,^]|. WEP tests can 
be conducted with laboratory or astronomical bodies be- 
cause the relevant differences are in the test-body com- 
positions. 

In its strong form the EP is extended to cover the 
gravitational properties resulting from gravitational en- 
ergy itself. In other words, it is an assumption about the 
way that gravity begets gravity, i.e. about the non-linear 
property of gravitation. Although general relativity as- 
sumes that the SEP is exact, alternate metric theories 
of gravity such as those involving scalar fields, and other 
extensions of gravity theory, typically violate the SEP 



[U0]. For the SEP case, the relevant test body differ- 
ences are in the fractional contributions to their masses 
by gravitational self-energy. To facilitate investigation of 
a possible violation of the SEP, the ratio between gravi- 
tational and inertial masses, Mq/Mj, is expressed as 
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where U is the body's gravitational self-energy (U < 0), 
Mc 2 is its total mass-energy, and 77 is a dimensionless 
constant [j| . Because gravitational self-energy is propor- 
tional to M 2 (i.e. U /Mc 2 oc M) and gravity is so ex- 
tremely weak, SEP test bodies that differ significantly 
must have astronomical sizes. Currently, the Earth- 
Moon-Sun system provides the best arena for testing the 
SEP with LLR being the only solar system technique 
available to enable the tests. 

The quasi-Newtonian acceleration of the Moon (to) 
with respect to the Earth (e), a = a m — a e , for the three- 
body Earth-Moon-Sun (s) system is 
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where fi* = fi e (M G /Mi) m + n m {M G /Mi) e and n k = 
GMk- The last two terms of Eq. J2J represent the solar 
effect on the motion of the Moon with respect to the 
Earth. A violation of the EP would produce a lunar 
orbit perturbation proportional to the difference in the 
two Mq /Mi ratios. An SEP test would involve Q 
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In general, rj is a linear function of seven of the ten PPN 
parameters, but considering only /3 and 7 



7] = 4/3 - 7 - 3. 



(4) 



In general relativity r\ — 0. A non-zero value for rj 
would produce a displacement of the lunar orbit about 
the Earth 0,01; a un it value would cause a 13 meter 
monthly range modulation pll |22| . 



III. NEW LLR TESTS OF RELATIVITY 

Each observation used in this analysis is a measured 
round-trip light time, here called a "range," between an 
observatory and a retroreflector. For data processing, the 
ranges represented by the returned photons are statisti- 
cally combined into a normal point; each normal point 
comprising from 3 up to ^100 photons. This paper's 
analysis of the LLR data from March 1970 to April 2004 
uses a total of 15,553 LLR normal points from McDonald 
and Haleakala Observatories and Observatoire de la Cote 
d'Azur. For the last 10 years of ranges the weighted rms 
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scatter after the fits is <~ 2 cm. This scatter is 0.5 x 10~ 10 
relative to the 385,000 km mean distance of the Moon. 

All fits of the lunar laser ranges involve a number of 
standard solution parameters for the Earth, Moon and 
lunar orbit (see |3J, |23j for the relativistic model used 
for JPL solutions). The ephemerides for the Moon and 
planets plus the rotation of the Moon are generated by 
a simultaneous numerical integration. Least-squares so- 
lutions require partial derivatives of range with respect 
to all solution parameters. Partial derivatives for the lu- 
nar orbit and rotation variations with respect to solution 
parameters are generated by numerical integration. 



A. Equivalence Principle Solution 

In essence, LLR tests of the EP compare the free-fall 
accelerations of the Earth and Moon toward the Sun. If 
the EP is violated, the lunar orbit will be displaced along 
the Earth-Sun line producing a range signature having a 
29.53 day synodic period [25| (different from the 

lunar orbit period of 27 days). Since the first LLR tests 
of the EP were published in 1976 0, the precision of 
the test has increased by two orders-of-magnitude 0, @ • 

The EP test is sensitive to the difference in Mq/Mi 
between the Earth and Moon. A test of the EP, corrected 
for solar radiation pressure p6| , is obtained from a fit of 
LLR data 
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This is equivalent to an orbit perturbation of Ar = 
(2.8 ± 4.1) mm cosD, where angle D corresponds to the 
29.53 day mean period of the new-full-new Moon cycle. 



B. Equivalence Principle Implications 

The LLR result © is a strong test of the EP. This test 
is sensitive to violations due to composition and gravi- 
tational self-energy. A University of Washington (UW) 
laboratory EP experiment [27I is designed to simulate the 
compositional differences of the Earth and Moon. That 
test of the relative acceleration is (1.0 ± 1.4) x 10~ 13 , 
where systematic and random uncertainties are combined 
[27| | . The laboratory results are insensitive to self-energy. 
A combination of the UW composition test with the LLR 
result (Eq. |SJ) yields the following result for the SEP test 
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= (-2.0±2.0) x 10~ 13 . (6) 
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Tests for violations of the EP due to self-energy are 
sensitive to a linear combination of the PPN quantities, 
Eq. Q. Considering only PPN (3 and 7, combine Eqs. © 
and 10 to obtain 



This expression would be null for general relativity, hence 
the small value is consistent with Einstein's theory. 

The SEP relates to the non-linearity of gravity (how 
gravity affects itself), with the PPN parameter (3 repre- 
senting the degree of non-linearity. Thus, LLR provides 
great sensitivity to (3, as suggested by the strong depen- 
dence of ij on (3 in Eqs. I@J and 10. The parameter 7 has 
been measured independently via time-delay and gravita- 
tional ray-bending techniques. The published Viking pj 
and Very Long-Baseline Interferometry (VLBI) Q uncer- 
tainties for 7 of ^0.002 led to a (3 uncertainty of ~0.0005 
dominated by the uncertainty in 7 [||. 

A much more accurate result for 7 was recently re- 
ported by the Cassini experiment; the test provided 
a verification that 7 is unity to a very high accuracy 
7 - 1 = (2.1 ± 2.3) x 10" 5 This leads to a significant 
improvement in the parameter f3 derived from r\. Com- 
bining the Cassini spacecraft determination of 7 with the 
i] of Eq. 10, determined from the LLR test of the EP and 
laboratory WEP results, gives 



13-1 = (1.2 ± 1.1) x 10" 



(8) 



This result is not a significant deviation of (3 from unity. 



LLR Tests of Other Gravitational Physics 
Parameters 



In addition to the SEP constraint Eq. 0, the PPN 
parameters 7 and [3 affect the orbits of relativistic point 
masses, and 7 also influences time delay Q. LLR can 
test this orbital f3 and 7 dependence, as well as geodetic 
de-Sitter precession, and G/G 0, S El. The possibil- 
ity of a time variation of the constant of gravitation, G, 
was first considered by Dirac in 1938 on the basis of his 
large number hypothesis, and later developed by Brans 
and Dicke in their theory of gravitation (for more details 
consult ^3 ) • Variation might be related to the expansion 
of the Universe, in which case G/G = (jHq, where Ho is 
the Hubble constant and a is a dimensionless parameter 
whose value depends on both the gravitational constant 
and the cosmological model considered. Revival of inter- 
est in the Brans-Dicke-like theories, with a variable G, 
is partially motivated by the appearance of superstring 
theories where G is a dynamical quantity .28'] . 

In this LLR analysis, the test of temporal variation of 
the gravitational constant results in 



G/G = (4 ±9) x 10~ 13 yr" 



(9) 



v = 4/3 - 7- 3= (4.4 ± 4.5) x 10~ 4 . 



(7) 



with a largest correlation of 0.74 with the diurnal tidal 
dissipation parameter. The G/G uncertainty is 83 
times smaller than the inverse age of the Universe, 
<o = 13.4 Gyr with the value for the Hubble constant 
Ho = 72 km/sec/Mpc from the WMAP data 29]. Any 
isotropic expansion of the Earth's orbit which conserves 
angular momentum will mimic the effect of G on the 
Earth's semimajor axis, a/a = —G/G 0. There is no 
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evidence for such local (~1 AU) scale expansion of the so- 
lar system. The uncertainty for G/G is improving rapidly 
because the sensitivity for new observations depends on 
the square of the time span. 

The test of geodetic precession yields 

K gp = -0.0019 ± 0.0064. (10) 

K gp is a relative deviation of geodetic precession from 
its GR value. The geodetic precession is highly corre- 
lated (0.88) with the lunar potential Love number and a 
parameter for lunar core oblateness. Adding the latter 
parameter, not present in the earlier solutions [ill 13. Il5|. 
increases the uncertainty of the geodetic precession. 

IV. CONCLUSIONS 



classical precession rates due to inclination and lunar Ji 
are much improved because of the added data span, but 
lunar core oblateness, a recent addition to the set of solu- 
tion parameters, is found to influence classical precession 
rates through its correlation with tidal distortion. The 
LLR EP test has improved markedly in the past decade, 
but much of that improvement was present in Q . Since 
that paper, the laboratory results for the WEP have im- 
proved by a factor of two and the determination of 
7 has improved by two orders-of-magnitude |10| yield- 
ing a five-fold improvement in PPN (3. Increased data 
span and future improved accuracy should continue to 
improve LLR tests of gravitational physics. The high ac- 
curacy LLR station being installed at Apache Point j^lj 
should provide major opportunities. 



The LLR data set provides sensitive tests of EP, PPN 
/?, geodetic precession, and G/G. There have been ma- 
jor improvements in the solution uncertainties since the 
1996 results Q. This improvement is partly due to an 
additional decade of high quality ranges and partly due 
to improvements in the model and data fits. The mat- 
ter of energy dissipation in the Moon, which previously 
was a limitation for G/G in is now much better un- 
derstood [2(|. For geodetic precession, the influences on 
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